Energy gaps in high- Tc superconductors: BCS after all? 
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A major impediment to solving the problem 
of high-Tc superconductivity is the ongoing con- 
fusion about the magnitude, structure and dop- 
ing dependence of the superconducting gap, Aq, 
and of the mysterious pseudogap found in under- 
doped samples^. The pseudogap opens around 
the (7r,0) antinodes below a temperature T* leav- 
ing Fermi arcs across the remnant Fermi surface^ 
on which the superconducting gap forms at T^. 
One thing that seems agreed is that the ratio 
2Ao/kBTc well exceeds the BCS value and grows 
with underdoping^'^, suggesting unconventional, 
non-BCS superconductivity. Here we re-examine 
data from many spectroscopies, especially Raman 
Big and i?2g scattering^'^, and reconcile them all 
within a two-gap scenario showing that the points 
of disagreement are an artefact of spectral-weight 
loss arising from the pseudogap. Crucially, we 
find that Ao(p), or more generally the order pa- 
rameter, now scales with the mean-field Tc value, 
adopting the weak-coupling BCS ratio across the 
entire phase diagram. 

We have long argued^ from thermodynamic and 
Knight-shift data that there are two distinct energy 
gaps in the electronic density of states (DOS). These 
are plotted in Fig. [Tfa) from the specific heat data for 
Yo.8Cao.2Ba2Cu307_^. (Similar results were found for 
Bi2Sr2CaCu208+^ (Bi-2212)). Aq is the projected max- 
imum in the T=0, d-wave SC gap, A(k), at the (7r,0) 
antinode. Aq(p) falls monotonically with increasing dop- 
ing, and vanishes at about p = 0.27 holes/Cu at the 
superconductor/metal transition. Eg is the magnitude 
of the pseudogap (PG) at the antinode, where, at the 
time, we assumed a triangular gap in the DOS which 
fills with increasing temperature. This filling effectively 
reflects the growth of the Fermi arcs with increasing 
temperature^ii^. Eg also decreases with doping but, un- 
like Ao, falls abruptly to zero at Pcrit = 0.19. Subsequent 
studies^ show that the Eg {p) doping dependence is prob- 
ably sublinear with Eg{p) cx (pcrit — p)^'^. 

It is important next to consider the k-dependence of 
these gaps. The SC gap has the d-wave form, A{6) — 
Ao cos(20) with nodes at 0=45° where 6 is the angle 
around the Fermi surface (FS). In contrast, as shown 
by specific heat^, Raman^, NQRi°, NMR", STM^^ and 
ARPES*^, the ground-state (GS) PG is not nodal but 
closes at 9q (< 45°) leaving remnant Fermi arcs around 
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where Oq defines the end of the Fermi arc. It is on the 
remnant Fermi arc that the SC gap opens. With decreas- 
ing doping the GS Fermi arc narrows and only closes at 
the node as p — > 0.05 where superconductivity disap- 
pears. Here there is no remnant arc on which the SC gap 
may appear and the loss of superconductivity at that 
point follows naturally. 

We plot in Fig.[2]the doping and angular dependence of 
A and Eg . The GS magnitudes are averaged from all the 
data shown in Fig. [TJ The figure also shows (green curve 
on the zero-energy plane) the evolution of Oq with doping 
as deduced previously^, ranging from 6*0 = at p = 0.19 
where the PG first opens, to 6*0 = 45° at p = 0.05. The 
black curve is the contour of equality of the two gaps. 
[Note that at p = 0.15 the two gaps are equal at the 
antinode but the PG has removed all the states below the 
dashed blue PG curve that would otherwise be pushed up 
above the SC gap. Thus there is a loss of spectral weight 
that is reflected e.g. in a marked diminution of the T = 
superfluid densityi^. This persists all the way out to 
p = 0.19, where the PG finally closesi^. Importantly, this 
observed loss of GS superfluid density illustrates that the 
PG is unrelated to thermal phase fluctuations because at 
T = all thermal fluctuations vanish.] 

The impact of this gap structure on the DOS is il- 
lustrated by the calculated DOS for Tl2Ba2Cu06+5 (see 
below) shown in Fig. [3l Three cases are shown: with the 
PG only, the SC gap only, and both coexisting. The loss 
of spectral weight in the coherence peak due to the PG 
is evident, as is the shift in the SC gap feature from Aq 

with no PG to 



E^ with the PG. The calculated 



DOS (red curve) corresponds in detail to the tunneling 
DOS observed in recent STM datai^. The van Hove sin- 
gularity (vHs) can be seen at £' — Ep « -85 meV. 

This is our starting point. The question we address 
here is: are other techniques consistent with this phe- 
nomenology? As a general principle the PG is best ob- 
served above Tc while the SC gap, in the absence of fluc- 
tuations, is observed only below Tc. If, and when, the PG 
and SC gap coexist below Tc care is required to distin- 
guish them. As a visual aid, in all data shown in Fig. [T] 
gaps measured above Tc (the PG) are shown in green, 
while gaps measured only below Tc are shown in red. 
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FIG. 1: The doping dependence of the maximum d-wave SC gap, Ao, and the pseudogap. Eg. Green data points show the gap 
observed above Tc (i.e. the pseudogap) while red data points show the gap observed below Tc (with care, the SC gap), (a) 
from specific heat measurements for Yo.8Cao.2Ba2Cu307_5; (b) from infrared c-axis conductivity for Ro.8Cao.2Ba2Cu307_4; 
(c) from Raman Big and B2g measurements on three 90K superconductors; (d) from intrinsic tunneling, (e) from SIS/SIN 
tunneling, and (f) from ARPES. The open circles in panel (d) are values of 2.5kBT^-^ . 



We turn then to our first comparison, that of the c- 
axis infrared conductivity. We show in Fig. [IJb) the two 
distinct gaps observed in cUipsometry studiesi^ii^. The 
green data points show the PG determined from the loss 
of spectral weight below a frequency 2Eg/h, which begins 
already around 300K, while the red data points show the 
SC gap determined from the loss of spectral weight below 
a frequency 2Ao/7i which begins only below Tc- In the 
heavily overdoped region only the SC gap is observed. 
Most notably there is a range of doping in which both 
gaps are clearly seen to coexist, and where the relative 
magnitudes of Eg{p) and Aq{p) swap over. Both gaps 
closely match the specific- heat-derived gaps in panel (a) . 
This is very encouraging, but when we turn to Raman 
data the result is rather difi'erent, at least superficially. 

The recent revival of the two-gap scenario for cuprate 
superconductivity is due largely to the insights derived 
from Raman scattering in the Big and B2g modes^i^. But 
while the Big gap was found to decrease monotonically 
from a large magnitude comparable with the exchange 
interaction, J, the B2g gap more or less followed the 
dome-shaped Tc(p) phase curve. Typical data is shown 
in Fig. for Hg-120l6, Bi-2212i6 and Y-123i^, each 
having similar maximum Tc values. Again, the Big gaps 
are shown in green and red depending on whether they 
are observed above or only below Tc, respectively. The 
B2g gap is shown in black diamonds. Like the infrared 
gaps these are pairbreaking gaps so we have divided by 
a factor 2. The following features are apparent: (i) the 



Big gap, in its p-dependence, is like Eg in panels (a) 
and (b) but is actually larger; and (ii) as noted, the i?2g 
gap is not monotonic but follows the SC dome seen in 
Tc{p). These features led the Sacuto group to interpret 
the underdoped Big gap as the PG and the i?2g gap as 
a direct measure of the SC gap^!^. Hiifner et alJ^, in an 
influential review, draw the same conclusion. Thus, at 
face value, the Big and B2g gaps appear to contradict 
the infrared and specific heat data. However, the analy- 
sis assumes a single-gap model to interpret spectra in a 
two-gap scenario. Below, we reconcile all these data. 

As shown in the inserts to Fig. |DJa) and (b), the Big 
scattering symmetry probes around the gap antinodes 
at (tTjO), while B2g scattering is dominated by contribu- 
tions near the d-wave gap nodes. Where there is a single 
rf-wave SC gap the Raman scattering peak in Big is a 
pairbreaking peak at 2Ao and at ^v^Aq for B2g. The 
\/2 factor is for a circular FS. For a typical cuprate FS 
it is w 4/3. On the other hand, where there is a second 
coexisting gap centred on the antinodes, as in the case of 
a PG with finite Fermi arcs, this removes the pile-up of 
states above the SC gap near (7r,0) that contributes to the 
2Ao pair-breaking peak, pushing the Big gap to higher 
energy and the B2g gap to lower energy^. This merely 
refiects the change in DOS shown in Fig. [S] The scatter- 
ing intensity also falls. In short, the i?2g pair-breaking 
gap is not the SC gap, Aq and it is somewhat fortuitous 
that it roughly follows the Tc{p) phase curve. 
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FIG. 2: Schematic diagram, based on Fig. [T] showing the 
doping- and angular-dependence of the combined gap struc- 
ture for the SC gap, Ao(S) (red curves), and the pseudo- 
gap, Eg{6) (blue curves), where 6 is the angle around the 
Fermi surface. The dotted blue curves show the ground-state 
pseudogap behavior in the absence of superconductivity. The 
green curve on the zero-energy plane is the contour of doping 
dependent ground-state Fermi arcs. The black curve shows 
the contour where Eg{6,p) — A(9,p). 
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FIG. 3: The calculated DOS based on a rigid ARPES-derived 
dispersion for Tl2Ba2Cu06 for three cases: with PG only 
(green), with SC gap only (black), and with both (red) where 
the DOS acquires two gap features and loses spectral weight 
in the coherence peak. The sub-gap feature marked A'q causes 
the -B29 gap and approximates the order parameter. The vHs 
lies at E - Ef ~ -85 meV. 



We have calculated the detailed Raman Big and B2g 
response and the results are summarised in Fig. [J] The 
imaginary part of the unscreened non-resonant Raman 
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where the integral is over occupied states below Ep, 
A(k) — i Ao (cos kx — cos ky) is the d-wave SC gap func- 
tion and E{\i) = ^e(k)2 -|- |A(k)|2. In the Big scattering 
symmetry 7(k)^i9 — jBig {coskx — cos ky). This func- 
tion is plotted in the inset to Fig. IUJa) and, as noted, is 
maximal on the antinodal sections of the FS. For B2g, 



7(k)^ 



'jB2g sin kx sin ky and the response, shown in 



the inset to Fig. [Jb), is mainly nodal. To compare 
with the HgBa2Cu05 data in Fig. [TJc), we use a rigid 
dispersion reported for Tl2Ba2Cu06 from fitted ARPES 
datcfc^S. The PG is modeled with doping-dependent Fermi 
arcs using Eq. O. We fit the parameter values for Eg (p) 
and 9o by assuming Aq and reproducing the Big and i?2g 
gaps, respectively, shown in Fig. (He). We also calculate 
spectra in the absence of the PG. The hole concentration 
is obtained from integrating the DOS. 

Fig. d] shows the calculated Big and i?2g spectra with 
(panels (a) and (b)) and without ((c) and (d)) the PG. 
Leaving aside the Raman continuum all features corre- 
spond in detail to those reported in Fig. [Ijc) and in the 
raw spectral^, including the reduction in peak intensity 
for B2g- In the absence of the PG both modes pro- 
gressively shift to higher energy with decreasing doping 
whereas, with the PG, the peak in the i?2g mode rises 
then falls (and weakens) as shown by the dashed arrow. 
This reflects the impact of the PG in removing antinodal 
spectral weight. At the same time the Big gap continues 
to rise, though more rapidly than in the absence of the 
PG. The peak energies are plotted in Fig UJe) and the 
fitting parameters Eg and 6*0 are plotted in Fig[4l[f). 

We note the following trends. In the absence of the 
PG the Big and B2g gap features track each other with 
a ratio of 4/3, as expected, and they precisely track Aq as 
determined from the peak-to-peak gap in the DOS. With 
the introduction of the PG the Big antinodal gap feature 
continues to take the magnitude of the SC gap in the 
heavily overdoped region beyond Pcrit , but below Pcrit it 
rises rapidly and its magnitude actually exceeds Eg just 



as we have noted. Here the Big gap is 2y Aq + -Eg which 

is plotted by the black symbols in Fig. Ufe) . 

The vHs is seen in panels (a) and (c) by the second 
(high-energy) peak in Big. Experimentally, this would 
be concealed by the electronic continuum above the gap 
but could emerge from the electronic background at high 
doping near or beyond the vHs. A possible observation 
of the vHs in Big Raman scattering at high doping is 
reported by Limonov et al.^^ . 

Turning to the B2g gap it is notable that this trends to 
zero at the SC/insulator boundary at p 0.05 and again 
one can understand how the B2g gap could be confused 
with Aq. But it falls to zero because, as shown, the 
ground-state Fermi arc shortens with decreasing doping, 
collapsing to the node at p « 0.05. The residual DOS 



4 




FIG. 4: The calculated Raman response for HgBa2Cu04 for various doping levels listed in the legends of panels (c) and (d): 
(a) and (b) with the PG, (c) and (d) without the PG. Insets show 7(k) from Eq.[2]plotted in the first quadrant of the Brillouin 
zone for Big (a) and for B2g (b). Solid black curves are the Fermi surface. Dashed arrows in (b) and (d) contrast the rise and 
fall in B2g gap in the presence of a PG with the monotonic behaviour without a PG. Panel (e) shows Big (red) and -829 (blue) 
pair-breaking gap values with (solid symbols) and without (open symbols) the PG. Arrows indicate the effect of the PG. Black 
symbols are 2^ A'g + Eg. Open black circles with error bars are values of i.^ksT^^ divided by 4/3, the renormalisation factor 
between the Big and -B29 gaps, (f) Parameter values obtained from the data fits. Crosses show 9o determined from STMi^. 



pile-up is then pushed all the way to the node, weakening 
as it approaches, and vanishes there. This result is borne 
out again in the present data fits - see Fig. [^f) where 
Oq 45° asp ^ 0.05. Note that if the GS of the PG were 
always nodal then our calculations show that the B2g 
gap follows the Big gap (reduced by the renomalisation 
factor of 3/4). It is the non- monotonic behavior of the 
gap that confirms the finite Fermi arc in the GS 
PG. The B2g pair-breaking gap can be approximated by 

I'^Aq — (|i?g)^ and roughly follows the magnitude of 

the SC order parameter, Aq. In general Aq < Aq but in 
the absence of the PG they are equal. 

The deduced values of Aq, Eg and Oq closely match 
those derived from specific heat. For 6*0, compare the 
green curve in Fig. [2] with blue diamonds in Fig. [H^f). 6*0 
values extracted from recent STM data on Fermi arcsiS, 
show the same trend (crosses in Fig. SJf)). We conclude 
that there is remarkable consistency between the specific 
heat, NMR Knight shift, c-axis infrared conductivity and 
Raman Big and i?2g scattering. 

We now turn to tunneling gaps - firstly intrinsic tun- 
neling. The results of intrinsic tunneling would appear 
to be regrettably downplayed. Typically involving up to 
20 bilayer stacks, intrinsic tunneling is arguably a bulk 
tunneling technique in comparison with scanning tunnel- 
ing spectroscopy which just probes the outer Cu02 layer. 
Intrinsic tunneling consistently reveals very sharp coher- 
ence peaks^^i^^ and the presence of two gaps with very 
different doping and temperature dependences. For ex- 
ample, investigations on Bi2(Sr2-2;La2;)Cu06+5 reveal a 



SC gap that closes at Tc with 2Ao/fcsTc 4.2 while the 
pseudogap remains fixed in value with increasing temper- 
ature and fills rather than closes^. With increasing dop- 
ing the pseudogap reduces in magnitude and falls to zero 
at p « 0.20 holes/Cu. This phenomenology is wholly con- 
sistent with what has been described above. In the case of 
Bi-2212 the evolution of Eg and Ao with doping has been 
reported in detail by Krasnov and coworkersSS. This sys- 
tem shows a similar behavior to Bi-2201 discussed above 
and the gap values are plotted in Fig. [Tfd) by the solid 
symbols. Here each gap can be discerned, both when 
Eg > Aq and when Eg < Aq. The detailed variation 
with doping is again consistent with all the data shown 
in panels (a) to (c), with the pseudogap falling abruptly 
to zero atp — 0.19. Such data has been questioned on the 
basis of overheating of the nanoscale mesas2i but this can 
be addressed and eliminated^'^'^*'. More recent intrinsic 
tunneling studies by Bcnscman et ali^ have been carried 
out in closely-spaced increments in doping. Their values 
of Aq are also plotted in Fig. [Hd) by the open symbols. 
They fully confirm the Krasnov data. 

SIN and SIS break-junction tunneling spectra have 
been measured for Bi-2212 by Miyakawa et alM-^ and 
their low-T gaps are plotted by the solid symbols in 
Fig. [H^e). These authors drew attention to the very 
large gaps observed in the underdoped region where 
2Ao/kBTc was found to progressively grow with under- 
doping (reaching 18.9 for a Tc=70K sample). Because 
these gaps tended to decrease near Tc they were pre- 
sumed to be SC gaps which profoundly exceeded BCS 
weak-coupling behavior. Unlike all other panels in Fig. [1] 
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here the green and red symbols do not represent gaps 
reported above and below Tj, respectively. No clear PG 
feature was observed above T^- However the data do 
show the same break in slope at p ~ 0.13 consistent with 
the crossover from PG to SG gap, and retaining the col- 
oration assists the eye to see this. These authors also 
noted that their gaps followed the doping dependence of 
the ARPES (7r,0) gaps. But we now know the ARPES 
(tTjO) gaps in the underdoped region to originate from 
the PG. In fact the conductance peaks show exactly the 
behavior expected from the square of the DOS shown 
in Fig. [31 The coherence peaks are large for the over- 
doped samples where the sole gap is the SC gap, Aq. But 
for the underdoped samples they have shrunk markedly 
(see Fig.l of rcf('*). The break occurs precisely at the 
crossover from red to green data points in Fig. [TJe). All 
four green data points have weak coherence peaks, as ex- 
pected f or a coexisting PG. M oreover, for the dispersion 
_Ek = A(k)2 + Eg{'ky, even if the pseudogap is 

the dominant gap, the composite gap feature will move 
to lower energy as T — > Tc because Ao(T) 0. Such a 
partial contraction of the gap energy near is not an 
indication the gap is a purely SG gap. 

Thus we feel that the general features seen in SIN and 
SIS tunneling are rather consistent with all the other 
spectroscopic results noted above. Moreover, the picture 
we present is entirely consistent with recent STM tunnel- 
ing datai^. What is puzzling is why the PG is not clearly 
seen above Tc in SIN and SIS. This could be a combina- 
tion of the k-dependent tunneling matrix elements and 
the fact that the (7r,0) quasiparticles, which are coherent 
in the SG state, become incoherent in the normal state. 
Below Tc, we have previously^ drawn attention to addi- 
tional second-gap features in the tunneling spectra that 
could be the pseudogap {Eg < Aq) in the case of the 
p = 0.16 sample'^, and the SG gap (Aq < Eg in the case 
of the p = 0.112 sample^. These are reminiscent of the 
second-gap features seen in Fig. [3] and are shown by the 
open green and red symbols, respectively, in Fig.^e). 

Our last comparison is with ARPES measurements of 
energy gaps. There is a very large body of ARPES litera- 
ture reporting energy gaps but this is complicated by the 
fact that many of these are mid-point, leading-edge gaps 
which tend to underestimate the gap and many of the re- 
mainder are peak-to-Ep energy differences which poten- 
tially overestimate the gap. The most reliable approach 
is to divide out the Fermi function, symmetrize the spec- 
tra and determine the peak-to-peak gap. Even then this 
does not necessarily differentiate between the PG and 
the SC gap. This requires either measurements above 
and below Tc or k-dependent measurements around the 
Fermi surface, or preferably both. Two papers from the 
Shen group^*'^^ for Bi-2212 do exactly this and thus en- 
able clear separation of the pseudogap and the SC gap. 
The projected (7r,0) gaps are plotted as a function of dop- 
ing in Fig. [TJf). While the scatter is quite large the same 
trend is seen as in all the other spectroscopies we have 
considered, with the PG overtaking the magnitude of the 



SG gap at about p « 0.13. We also plot the data of Gam- 
puzano et al.^^ which is the simple low-T quasiparticle 
peak position. It therefore represents the combination of 
both gaps when they are both present. It also shows the 
same break or kink in its doping evolution. 

Finally, we compare the Aq values with the mean-field 
Tc value, T^f, determined fr om the fluctuation specific 
heat. We have shown^, that both amplitude and phase 
fluctuations set in simultaneously above the observed Tc 
thus shifting Tc well below T™-^. Using an entropy- 
balance procedure we found that even in the overdoped 
region this downward shift exceeds lOK. But in the un- 
derdoped region the shift rises to as much as 60 or 70K 
in the case of Bi-2212 (see Tc data in Fig. [DJd)). Sig- 
nificantly, while 2A/kBTc is variable and grows to large 
values at lower doping, 2A/kBT™f remains constant and 
comparable to the d-wave weak-coupling BCS value, 4.3. 
It is T"''f which is the true SG energy scale, not Tc- In 
Fig. [TJd) we plot values of 2.3 x T^f and they track the 
Aq values until Eg becomes comparable to Aq when the 
pseudogap is large enough to significantly reduce the or- 
der parameter and T™^ then falls. Indeed, as noted, the 
B2g gap is an approximate measure of (2 x |x) the SG 
order parameter, Aq (not, as stated, to be confused with 
Aq unless Eg = 0). As shown by the open black circles 
in Fig. mje) the quantity 4.3xTc™^ follows the B2g gap 
rather closely, implying that 2AQ/fcsT™^ remains close 
to the weak-coupling BCS value 4.3 across the entire dop- 
ing range. This is a central result. A similar comparison 
with the specific heat or infrared Aq values in Fig. [TJa) 
or (b) also returns values of 2Ao/kBT™f very close to 
4.3 when Eg is small or absent. 

The presence of SC fluctuations will cause a partial 
tunneling gap above T^ which must not be confused 
with the PG. It too will have washed out coherence peaks 
but will only persist to T™-'^ which is usually much less 
than the PG temperature, T*— . It is SC fluctuations that 
give the impression that the SC gap evolves smoothly 
into the PG'^^'^'^ and it is the SC fluctuation gap extend- 
ing only lOK or so above Tc that has been observed in 
overdoped STM data^^. In both cases this has led to con- 
fusion between the SG gap and the PG. It is the fact that 
Aq and Eg, and Tc, T^^ and T* are all comparable in 
magnitude^ , combined with the presence of fluctuations 
above Tc that has sustained this long-standing confusion 
which hopefully can now be laid to rest. 

In conclusion, we have examined specific heat, c-axis 
infrared spectroscopy, Raman Big and B2g scattering, in- 
trinsic tunneling, SIS/SIN & STM tunnefing and ARPES 
in HTS cuprates. We find that they all present a consis- 
tent picture of two monotonically increasing energy gaps 
as doping is reduced. The pseudogap. Eg, rises rapidly 
from zero at Pcrit — 0.19 and reaches a scale of ~ J at 
p = 0, where J is the exchange energy. It forms first 
around (7r,0) leaving ungapped Fermi arcs on which the 
SG gap opens. While the Fermi arcs decrease with de- 
creasing temperature they still remain finite in length 
at T = 0. The arcs also shorten with decreasing dop- 
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ing, pinching off at the d-wave nodes at p « 0.05 where 
SC disappears. The T = SC gap opens at p « 0.27 
and rises less rapidly but always monotonically, reach- 
ing a scale of about ^J. The B2g Raman gap, which 
naively shows an apparent scaling with the dome-shaped 
Tc curve must be reinterpreted in a two-gap scenario and 
the downturn is shown to be due to redistributed spec- 
tral weight. When corrected, both Big and B2g Raman 
scattering present a consistent picture of a monotonically 
increasing SC gap. Unfortunately, Hiifner et alf^ reiter- 
ate this naive picture of a dome-shaped B2g SC gap in 
their review, paying no attention to the long-standing 



thermodynamic, NMR and (more recent) infrared data. 
Their resultant gap plot is further reproduced in detail 
by Kohsaka et al.^'^ and has quickly become widely ac- 
cepted. We hope the present work refutes this picture. 
Most importantly, Aq does not scale with Tc but with its 
mean-field value, T™^ across the overdoped region, with 
2Ao/A:bTJ"-^ remaining close to the BCS weak-coupling 
value of 4.3. With the opening of the pseudogap the 
order parameter Ag falls below Aq but it seems that, 
quite generally, 2Aq/ ksT^^ « 4.3 at all doing levels. 
The cuprates may be more conventional than we once 
thought. 
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